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HI power spectrum of the spiral galaxy NGC628 
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ABSTRACT 

We have measured the HI power spectrum of the nearly face-on spiral galaxy NGC628 
(M74) using a visibility based estimator. The power spectrum is well fitted by a power 
law P(U) = AU a , with a = - 1.6 ± 0.2 over the length scale 800pcto8kpc. The 
slope is found to be independent of the width of the velocity channel. This value of the 
slope is a little more than one in excess of what has been seen at considerably smaller 
length scales in the Milky- Way, Small Magellanic Cloud (LMC), Large Magellanic 
Cloud (SMC) and the dwarf galaxy DDO210. We interpret this difference as indicating 
a transition from three dimensional turbulence at small scales to two dimensional 
turbulence in the plane of the galaxy's disk at length scales larger than galaxy's HI 
scale height. The slope measured here is similar to that found at large scales in the 
LMC. Our analysis also places an upper limit to the galaxy's scale height at 800 pc . 

Key words: physical data and process: turbulence-galaxy:disc-galaxies:ISM 
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1 INTRODUCTION 

Evidence has been mounting in recent years that turbulence 
plays an important role in the physics of the ISM as well 
as in governing star formation. It is believed that turbu- 
lence is responsible for generating the hierarchy of struc- 
ture s present across a range of spatial scales in the ISM 
fe.g. lElmegren fc scald kuWal ; lElmegren fc scald l2004bl ). In 
such models the ISM has a fractal structure and the power 
spectrum of density fluctuations is a power law, indicating 
that there is no preferred "cloud" size. 

On the observational front, the power spectrum anal- 
ysis of HI intensity fluctuations is an important tech- 
nique to probe the structure of the neutral ISM in galax- 
ies. The power spectrum analysis in our own Galaxy, 
Large Magellanic Cloud(LMC) and Small Magellanic 
Cloud(SMC) finds the power sp ectrum of the HI inten 



sity fluctuat i on to be a power lawdCrovisier fc Dicke 



Greenlll993l; Dcshpande et al. 2000; lElmegreen et all 
Stanimirovic et al. 1999) which is a characteristic of a turbu 
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lent medium. Similarly, Westpfahl et al.(1999) showed that 
the HI distribution in several galaxies in the M81 group is a 
fractal. Further, Willett et al.(2005) used the Fourier trans- 
form power spectra of the V and Ha images of a sample of 
irregular galaxies to show that the power spectra in opti- 
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cal and Ha pass-bands are also power law, indicating that 
there is no characteristic mass or luminosity scale for OB 
associations and star complexes. 

Recently iBegum et alj l|2006l) presented a visibility 
based formalism for determining the power spectrum of HI 
intensity fluctuations in galaxies whose emission is extremely 
weak. This was applied to a dwarf galaxy, DDO 210. In- 
terestingly, the HI power spectrum of this extremely faint, 
largely quiescent galaxy was found to be a power law with 
the same slope as that observed in much brighter galaxies. 
In this paper the same formalism is used to measure the 
power spectrum of HI intensity fluctuation in the nearby 
spiral galaxy NGC628. 

NGC628(M74) is a nearly face-on SA(s)c spiral 
galaxy with an inclinatio n angle in the range 6° to 13° 
ijKamphuis fc Briggsl fl993 >. It has a very large HI disk ex- 
tending out to more than 3 times the Holmberg diameter 
|Kamphuis fc Briggsl [l992T) . Elmegreen et al. (2006) have 
found a scale-free size and luminosity distribution of star 
forming regions in this galaxy, indicating turbulence to be 
functional here. The distance to this galaxy is uncertain 
with previous est imat es ranging from 6.5 Mpc t o 10 Mpc. 
iBriggs et al] Il980l and iKamphuis fc Briggsl Il992l have used 
a Hubble flow dista nce of ~10 Mpc. On the other hand, 
ISharina et alj (|l996l ) estimated a distance of 7.8 ± 0.9 Mpc 
from the brightest blue star in the galaxy. This distance 
estimate matche s with an independent photometric dis- 
tance estimate bylSohn fc Davidge I ll 19961 ) . In a recent study 
IVinko et all l|2004h inferred the distance to be 6.7±4.5 Mpc 
by applying the expanding photosphere method to the hy- 



2 Prasun Dutta, Ayesha Begum, Somnath Bharadwaj and Jayaram N. Chengalur 



6 8 10 




01 34 45 30 15 00 33 45 30 15 

RIGHT ASCENSION (B1950) 



Figure 1. The 38" X 36" resolution integrated HI column density 
map of NGC628 (contours) overlayed on the optical DSS image 
(grey scale). The contours levels are 0.24, 1.96, 3.68, 5.41, 7.13, 
8.85, 10.57 and 12.29 xlO 20 cm" 2 . 

per novae SN2002ap. Throughout this paper we adopt the 
photometric distance of 8 Mpc for NGC628. At this distance 
1" corresponds to 38.8 pc. 

In this paper we present the power spectrum of HI in- 
tensity fluctuations of NGC62 8 derived using the visibility 
based formalism developed by IBegum" et al.l (|2006l ). Studies 
of external galaxies like NGC628 which has a very extended 
HI disk holds the potential of probing the ISM and its power 
spectrum at leng th scales m uch larger than has been possible 
in earlier studies ( Oev , 2 0021 ) . 



2 DATA AND ANALYSIS 

We have used archival HI data of NGC628 from the Very 
Large Array (VLA). The observations had been carried out 
on 1 st August and 14 t/l November, 1993 respectively in 
the C and D configurations of the VLA, as a part of the 
A AO 163 observing program. The multi-configuration data 
were downloaded from the VLA archive and reduced in the 
usual way using standard tasks in classic AIPS3. For each 
VLA configuration, bad visibility points were edited out, af- 
ter which the data were calibrated. The calibrated data for 
both configurations was combined using DBCON. The HI 
emission from NGC628 spans 64 central channels of the 256 
channel spectral cube (with channel width 1.29 km s _1 ). 
A continuum image was made using the average of all the 
line free channels. The continuum from the galaxy was sub- 
tracted from the data in the uv plane using the AIPS task 
UVSUB. The resulting continuum subtracted data was used 
for the subsequent analysis. Figure Q] shows a total HI col- 
umn density (Moment 0) map of NGC628 from an image 

1 NRAO Astrophysical Image Processing System, a commonly 
used software for radio data processing. 



made from this data. The HI disc of the galaxy is nearly 
face-on. The angular extent of the HI distribution in Fig- 
ure [T] is roughly 1 1' x 15'. Using deep VLA D array mo- 
saic observations, dKamphuis fe Briggslll99"2h had detected 
a ~ 38' x 31' faint diffuse HI envelope around NGC628. This 
extended envelope is not detected in the current dataset; the 
emission that we do detect is instead restricted to the main 
HI disc of NGC 628. The results we discuss below are hence 
also relevant only to the gas in the main HI disc around 
NG C628. 

IBegum et alj l|2006l ) contains a detailed discussion 
of the visibility based HI power spectrum estimator 
Pm(U) =< V„(U)V*(& + AU) >, we present only a brief 
discussion here. Every visibility V V (U) is correlated with all 
other visibilities V V (U + AU) within a disk | AU \< D. We 
shall discuss the considerations for deciding the value of D 
shortly. The correlations are averaged over different U direc- 
tions assuming that the signal is isotropic. To increase the 
signal to noise ratio we further average the correlations in 
bins of U and over all frequency channels with HI emission. 
The expectation value of the estimator Pm{U) is real and 
it is the convolution of the HI power spectrum Pm{U) and 
a function |W/„([/)| 2 which can be assumed to be sharply 
peaked around U — with a width of order D. At baselines 
U S> D the function IW^t/)! 2 can be well approximated by 
a Dirac delta function and the expectation value of Phi (U) 
gives an estimate of the HI power spectrum Phi (U). The 
value of P> -1 is of the order of the angular extent of the HI 
emission (Figure[l|, and we use D = QAKX for our analysis. 

The 1 — a error-bars for the estimated power spectrum 
is a sum, in quadrature, of contributions from two sources of 
uncertainty. At small U the uncertainty is dominated by the 
fact that we have a finite and limited number of independent 
estimates of the true power spectrum, while at large U it is 
dominated by the system noise in each visibility. 

The HI emission spans 64 frequency channels each of 
width 1.29 km s _1 . To determine if the slope of the HI power 
spectrum changes with the width of the frequency channel, 
we have combined N successive channels to obtain a data 
set with 64/iV channels of width N x 1.29 kms -1 each. We 
have determined the HI power spectrum for a range of N 
values. 



3 RESULTS AND DISCUSSION 

Figure [2j\ shows the real and imaginary parts of Phi (U), 
which is the observed value of the estimator Pm(U) for the 
64 channels which have HI emission. As expected from the 
theoretical considerations mentioned earlier, the imaginary 
part is well suppressed compared to the real part. To test 
for a possible contribution from residual continuum, we also 
show the real part of Pm(U) using 64 line free channels. This 
is found to be much smaller than the signal. For the chan- 
nels with HI emission the observed Phi (U) may be directly 
interpreted as the HI power spectrum at U values that are 
considerably larger than D — 0.4 KA. We find that a power 
law Pm(U) — AU a with slope a — —1.6 ± 0.2 provides a 
good fit to the results over the U range 1.0 K A to 10.0 K A 
(Figure 03) which corresponds to spatial scales of 800pc to 
8kpc. 

Both HI density fluctuations as well as spatial fluctua- 
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Figure 2. A) Real and imaginary parts of the observed value of the HI power spectrum estimator Pm(U) for TV = 64 ie. all frequency 
channels with HI emission were collapsed into a single channel. The real part is also shown using 64 line-free channels collapsed into a 
single channel. B) Best fit power law to Pui{U). The channel width is varied (TV = 1,32,64; top to bottom), la error-bars are shown 
only for TV = 64. 



N 


AV 


a 


la 




(kms" 1 ) 
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1.3 


-1.7 


±0.2 


32 


42.3 


-1.6 


±0.2 


64 


82.5 


-1.6 


±0.2 



Table 1. NGC628 power spectrum are well fitted by a power law 
of Phi = AU a . This table sunrises the result for different channel 
width. 



tions in the velocity of the HI gas contribute to fluctuations 
in the HI specific int e nsity. Considering a turbulent ISM, 
lLazarian fc Pogosvanl (2000) have shown that it is possible 
to disentangle these two contributions by studying the be- 
havior of the HI power spectrum as the thickness of the fre- 
quency channel is varied. If the observed HI power spectrum 
is due to the gas velocities, the slope of the power spectrum 
is predicted to change as the frequency channel thickness is 
increased. We have tested this by determining the HI power 
spectrum for different values of the channel width in the 
range 1.29 kms" 1 to 82.6 kms -1 (Figure[5j3) and do not find 
any change in the slope of the HI power spectrum (Table [IJ. 
As the thickest channel that we have used is considerably 
wider than the typical HI velocity dispersion of 7 — 10 km s~ 
seen in spiral galaxies l|Shostak fc van der Kruitlll984T l. we 
conclude that the observed HI power spectrum of NGC628 
is purely du e to density fluctua tions. Our finding is simi- 
lar to that of I Begum et al.l (|2006t ) who noticed no change of 
the slop e with channel width for the dwarf galaxy DDO210. 
Further, lElmegreen et all l|200lf) also reported a similar be- 
havior for LMC. 

Earlier studies of the Milky- Way, and also of the dwarf 
galaxies LMC , SMC and DDO210 JCrovisier fc Dickevll983l; 



Greenl 1 19931; IStanimirovic et al.l 1 19991: iDeshpande et al.l 



200d : lElmegreen et all 1200 ll ; iBegum et all I2006T I have all 
found a power law HI power spectrum with slope ~ —3. 



On the contrary, we find a slope —1.6 ± 0.2 for NGC628. 
This is a little more than one in excess of the earlier 
values. However, when comparing these values it should 
be noted that the earlier works have all measured the 
HI power spectrum at much smaller length scales in the 
range 10 to 500 pc [MW (5pc - 200 pc), SMC(4pc - 
30 pc), DDO210(80pc — 500 pc)] whereas the current mea- 
surement probes much larger length scales from 800 pc to 
8.0 kpc. The typical HI scale heights within t he Milky- Way 
ijLockman et al. 1984; Wouterloot et al.lll990h and external 
galaxies (e.g. iNaravan fc Jodl2002l ) are well within 1.5 kpc 
. This implies that on the largest length scales which we 
have probed, the turbulence is definitely confined to the 
plane of the galaxy's d isk and is therefore two dimensional. 
lElmegreen et al] (|200lh have found that the HI power spec- 
trum of LMC flattens at large length scales, which was in- 
terpreted as a transition from three dimensional to two di- 
mensional turbulence. We conclude that the slope is differ- 
ent in our observations because it probes two dimensional 
turbulence, whereas the earlier observations were on length 
scales smaller than the scale height where we can expect 
three dimensional turbulence. To the best of our knowl- 
edge our results are the first observational determination 
of the HI power spectrum of an external spiral galaxy at 
such large length scales which are comparable to the radius 
of the galaxy's disk. 

IWestpfahl et al l (|l999h have performed a fractal analy- 
sis using the perimeter-area dimension of intensity contours 
in HI images of several galaxies in the M81 group. Of par- 
ticular interest is the galaxy M81, a spiral galaxy for which 
the perimeter-area dimension was found to be ~ 1.5 at a 
length scale ~ 10 kpc. This observation is consistent with a 
power law power spectrum of slope —1.6 ± 0.2 provided the 
assumption that the local dimension has the same value as 
the perimeter-area dimension is valid. 

It is difficul t to probe the HI sc ale height of external 
face-on galaxies. IPadoan et al.l (|200ll ) present a method to 
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probe the scale height from a chan ge in the slope of the Spe c- 
tral Correlation Function(SCF) [Rosolowskv et all l| 19991 )]. 
and applied it to HI dat a for the LMC to e stima te the scale 
height to be ~ 180pc. E lmegreen et all (|200ll ) suggested 
that one could use a change in the slope of the power spec- 
trum of the density fluctuations to measure the scale height 
of face on gas disks. Applying this method to HI data for 
the LMC they measure scale height of 100 pc. To the best of 
our knowledge, prior to this work, there has been no obser- 
vational constraint on the HI scale height of NGC628. Since 
the scale height is definitely less than 8kpc and the power 
spectrum is found to have the same slope from 800 pc to 
8kpc, we conclude that th e scale height must be less than 
800 pc. lKregel et al.l l|2004 ) present HI images of a large sam- 
ple of edge on intermediate to late type spirals; from their 
data the ratio of the HI disk height to the radius of the HI 
disk (at a column density of IMqpc— 2) is ~ 0.06 ± 0.15. 
From Fig. [Q the disk of NGC628 has a diameter of ~ 28kpc 
at this column density. From the average thickness to ra- 
dius ratio for edge on galaxies, one would expect NGC628 
to have a scale height of ~ 840pc, consistent with our re- 
sult. Future observations of NGC628 with higher angular 
resolution should be able to put a tighter constraint or even 
determine the scale height. 
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